Conservation laws, integrability and transport in one-dimensional quantum systems 
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In integrable one-dimensional quantum systems an infinite set of local conserved quantities exists 
which can prevent a current from decaying completely. For cases like the spin current in the XXZ 
model at zero magnetic field or the charge current in the attractive Hubbard model at half filling, 
however, the current operator does not have overlap with any of the local conserved quantities. We 
show that in these situations transport at finite temperatures is dominated by a diffusive contribution 
with the Drude weight being either small or even zero. For the XXZ model we discuss in detail the 
relation between our results, the phenomenological theory of spin diffusion, and measurements of the 
spin-lattice relaxation rate in spin chain compounds. Furthermore, we study the Haldane-Shastry 
model where the current operator is also orthogonal to the set of conserved quantities associated 
with integrability but becomes itself conserved in the thermodynamic limit. 



I. INTRODUCTION 

In classical dynamics the KAM theorem quantita- 
tively explains what level of perturbation can be exerted 
on an integrable system so that quasi-periodic motion 
survives. 1 A classical system with Hamiltonian H and 
phase space dimension 2N is integrable if N constants 
of motion Qk exist (i.e., the Poisson bracket vanishes, 
{H, Qk} — 0) which are pairwise different, {Qk, Qi} = 0. 
Defining integrability for a quantum system is, however, 
much more complicated and no analogue of the KAM 
theorem is known. For any quantum system in the ther- 
modynamic limit an infinite set of operators exist which 
commute with the Hamilton operator. This can be seen 
by considering, for example, the projection operators 
onto the eigenstates of the system, [H, \n)(n\] = with 
H\n) — E n \n). In quantum systems described by tight- 
binding models with short-range interactions it is there- 
fore important to distinguish between local conserved 
quantities Q n = ^2jq n ,j, where q n j is a density oper- 
ator acting on n adjacent sites j, and nonlocal conserved 
quantities like the projection operators mentioned above. 
In a field theory, a conserved operator is local if it can 
be written as integral of a fully local operator. Most 
commonly, quantum systems are called integrable if an 
infinite set of local conserved quantities exists which are 
pairwise different. This definition includes, in particu- 
lar, all Bethe ansatz integrable one-dimensional quantum 
systems. Here the local conserved quantities can be ex- 
plicitly obtained by taking consecutive derivatives of the 
logarithm of the appropriate quantum transfer matrix 
with respect to the spectral parameter. 2 

In recent years, many studies have been devoted 
to the question if integrability can stop a sys- 
tem from thermalizing 3-5 or a current from decaying 
completely. 6-19 That conservation laws and the trans- 
port properties of the considered system are intimately 
connected is obvious in linear response theory (Kubo 
formula), which relates the optical conductivity to the 



retarded equilibrium current-current Green's function. 
Specializing to a lattice model with nearest-neighbor hop- 
ping, the Kubo formula reads 76 
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Here, L is the system size, T the temperature, Skin the 
kinetic energy operator, J the spatial integral of the cur- 
rent density operator, and the brackets denote thermal 
average. The real part of the optical conductivity can be 
written as 



cr'((xi) = 2ttDS(uj) + Q- reg (uj). 



(1.2) 



A nonzero Drude weight D implies an infinite dc con- 
ductivity. Castella et al. [6] showed that it is possible to 
circumvent the direct calculation of the current-current 
correlation function in the Kubo formula and compute 
the finite temperature Drude weight using a generaliza- 
tion of Kohn's formula. 20 This formula relates the Drude 
weight to the curvature of the energy levels with respect 
to magnetic flux. The observation that integrable and 
nonintegrable models obey different level statistics, 21-23 
as well as the calculation of the Drude weight from exact 
diagonalization for finite size systems, 24 led to the con- 
jecture that anomalous transport, in the form of a finite 
D(T > 0), is a generic property of integrable models. 6 

This conjecture is corroborated by the relation between 
the Drude weight and the long-time asymptotic behavior 
of the current-current correlation function 7 
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In view of this relation, ballistic transport, D =/= 0, means 
that the current-current correlation function does not 
completely decay in time. The Qk operators in Eq. (1.3) 
form a set of commuting conserved quantities which are 
orthogonal in the sense that (QkQi) = {Q 2 k )o~ki- That 
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conserved quantities provide a lower bound for the long- 
time asymptotic value of correlation functions is a general 
result due to Mazur. 25 In fact, it can be shown that the 
equality holds if the right-hand side of Eq. (1.3) includes 
all conserved quantities Qk, local and non-local. 26 

The implications of Mazur's inequality for transport in 
integrable quantum systems were pointed out by Zotos 
et al. [7]. For integrable models where at least one con- 
served quantity Q n has nonzero overlap with the current 
operator, (Q n J) ^ 0, Mazur's inequality implies that 
the Drude weight is finite at finite temperatures. This 
happens, for instance, for charge transport in the Hub- 
bard model away from half-filling and for spin transport 
in the S = 1/2 XXZ model at finite magnetic field. 
By contrast, for nonintegrable models, for which all local 
nontrivial conservation laws are expected to be broken, 
the right-hand side of Eq. (1.3) is expected to vanish so 
that D(T > 0) = 0. In these cases the delta function in 
Eq. (1.2) is generically presumed to be broadened into 
a Lorentzian Drude peak. The dc conductivity is finite 
and the system is said to exhibit diffusive transport. 

However, in many cases of interest the known local 
conserved quantities associated with integrability are or- 
thogonal to the current operator because of their symme- 
try properties. This happens, for instance in the XXZ 
model at zero magnetic field h. In terms of spin- 1/2 op- 
erators S; the model reads 

N 

H = £[J(Sf Sf +1 + SfSf +1 + AS? S? +1 ) - hSf]. (1.4) 

Here N is the number of sites, A parametrizes an ex- 
change anisotropy, and J is the exchange constant which 
we set equal to 1 in the following. One can show that 
all local conserved quantities of the XXZ model are 
even under the transformation 5| — > —Sj, 5- — > Sj, 
whereas the current operator is odd. 7 In these cases the 
integrability-transport connection has remained a conjec- 
ture. Nonetheless, several works have presented support 
for a finite Drude weight in models where the current op- 
erator has no overlap with the local conserved quantities, 
at least in some parameter regimes. The list of meth- 
ods employed include exact diagonalization (ED), 14-16,27 
Quantum Monte Carlo (QMC) 12,28-30 and Bethe ansatz 

(BA) 8,9,31 

In particular for the XXZ model, ED results for chains 
of lengths up to L = 18 sites 16 suggest that at high 
temperatures the Drude weight extrapolates to a finite 
value in the thermodynamic limit for values of exchange 
anisotropy in the critical regime, including the Heisen- 
berg point. The Drude weight appears to vanish for large 
values of anisotropy in the gapped regime, but the min- 
imum value of anisotropy for which it vanishes cannot 
be determined precisely. Although no conclusions can be 
drawn about the low temperature regime, the claim is 
that if the Drude weight is finite at high temperatures 
it must also be finite and presumably even larger at low 
temperatures. The weakness of this method is that it 



assumes that the finite size scaling of D(T, L), which is 
not known analytically and is only obtained numerically 
for L < 18, can be extrapolated to the thermodynamic 
limit. This is not necessarily true for strongly interacting 
models where even at high temperatures there may be a 
large length scale above which the behavior of dynamical 
properties changes qualitatively. 

In the QMC approach for the XXZ model in Refs. [11, 
12], the Drude weight was obtained by analytic contin- 
uation of the conductivity function a(q,icu n ), which is 
a function of the Matsubara frequencies ui n , to real fre- 
quencies. This method uses a fitting function to try to 
extract a decay rate 7 which broadens the Drude peak 
if the Drude weight vanishes, but it clearly fails to find 
decay rates that are smaller than the separation between 
Matsubara frequencies, 7 <C T. In fact, the application of 
this method has even led to the conclusion that the Drude 
weight is finite for some gapless systems that are not 
integrable. 28 ' 29 This is hard to believe, considering that 
Eq. (1.3) requires the existence of nontrivial conservation 
laws which have a finite overlap with the current operator 
in the thermodynamic limit in order for the Drude weight 
to be finite. An attempt was made to explain the Drude 
weight for nonintegrable models described by a Luttinger 
liquid fixed point based on conformal field theory, 13 but 
this analysis neglects irrelevant interactions that lead to 
current decay and render the conductivity finite. 10 

Although the Drude weight for the XXZ model has 
been calculated exactly by BA at T = using Kohn's 
formula, 32 the calculation of D(T > 0) by BA is hindered 
by the need to resort to approximations in the treatment 
of the excited states. The BA calculation by Zotos 8 fol- 
lows an ansatz proposed by Fujimoto and Kawakami 31 
that employs the thermodynamic Bethe ansatz (TBA), 
which relies on the string hypothesis for bound states of 
magnons. 2 This approach predicts that the Drude weight 
is finite and decreases monotonically with T for the XXZ 
model in the critical regime at zero magnetic field, except 
at the Heisenberg point, where D(T) vanishes for all fi- 
nite temperatures. Benz et al. [9] criticized the TBA 
result and pointed out that it violates exact relations for 
D{T) at high temperatures. These authors presented an 
alternative BA calculation of the Drude weight based on 
the spinon and anti-spinon particle basis and predicted 
a different temperature dependence than the TBA re- 
sult. In particular, the Drude weight is found to be finite 
for the Heisenberg model at zero magnetic field. Ac- 
tually, for values of anisotropy near the isotropic point 
this approach predicts that D{T) increases with T at 
low temperatures. Like the TBA result, the result based 
on spinons and anti-spinons violates exact relations at 
high temperatures. A consistent calculation of the Drude 
weight by applying the BA to the finite temperature 
Kohn formula is therefore still an unresolved issue. 

Integer-spin Heisenberg lattice models are not inte- 
grable, but their low energy properties are often stud- 
ied in the framework of the continuum 0(3) nonlinear 
sigma model, which is integrable. Interestingly, BA cal- 
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culations for the nonlinear sigma model predict a finite 
D(T > 0) that is exponentially small at temperatures 
below the energy gap. 33,34 In fact, it has been argued 34 
that the finite Drude weight is due to at least one non- 
local conserved quantity of the quantum model which is 
known explicitly 35 and has overlap with the current oper- 
ator. This result is not consistent with the semiclassical 
results by Damle and Sachdev, 36 ' 37 which predict diffu- 
sive behavior for gapped spin chains at low temperatures. 
No nonlocal conserved quantities that overlap with the 
current operator are known explicitly for the integrable 
5 = 1/2 XX Z model. For finite chains such quantities 
can be constructed explicitly, however, in a numerical 
study of chains with L < 18 no definite conclusions could 
be made whether or not the overlap remains finite in the 
thermodynamic limit. 14 

The debate about the role of integrability in trans- 
port properties of integrable models is connected with 
the question of diffusion in S = 1/2 spin chains. The 
term spin diffusion first appeared in the context of the 
phenomenological theory, 38-40 where it refers to a char- 
acteristic form of long-time decay of the spin correla- 
tion function. In the phenomenological theory, the case 
where the total magnetization in the direction of quan- 
tization, S z = is conserved is considered. It 
is then said that spin diffusion occurs if the Fourier 
transform of the spin-spin correlation function G(q, t) = 
^ e-^-^- r °)(SJ(t)S^(0)) for small wavevector q de- 
cays with time as G(q, t) ~ e -£)q *, where D is the dif- 
fusion constant. Provided that the behavior at small q 
dominates, this implies that the Fourier transform decays 
as G(r,t) ~ i -d / 2 exp(-|r| 2 /4Di). This theory was for- 
mulated to explain results of inelastic neutron scattering 
experiments in three-dimensional ferromagnets at high 
temperatures. The assumptions of the theory are usually 
motivated by the picture that at high temperatures the 
spin modes are described by independent Gaussian fluc- 
tuations. The decay of the correlation function G(r,t) 
can then be interpreted as a random walk of the magne- 
tization through the lattice as in classical diffusion. 

It is important to note that this definition of diffu- 
sion is not obviously related to that of diffusive trans- 
port given earlier, since the two definitions refer to two 
different correlation functions. In the phenomenological 
theory, transport is found to be diffusive because the lo- 
cal magnetization obeys the diffusion equation and the dc 
conductivity is therefore finite. However, more generally 
diffusion in the autocorrelation function for the density 
of the conserved quantity does not exclude the possibility 
of ballistic transport, understood as a nonzero long-time 
value for the current-current correlation function. 

The applicability of the phenomenological theory of 
diffusion to spin chains described by the integrable XX Z 
model is of course questionable. Even at high tempera- 
tures, the spin dynamics is likely to be constrained by 
the nontrivial conservation laws and the assumption of 
independent modes is not expected to hold. In the one 
case where the long-time behavior of the autocorrelation 



function can be calculated exactly, namely the XX model, 
which is equivalent to free spinless fermions, diffusion 
does not occur since G(x = 0, t) <~ t -1 for large t, as 
opposed to i -1 / 2 expected for diffusion in one dimen- 
sion. For decades, a great deal of effort has been made 
to compute the autocorrelation function for the XXZ 
model with general values of anisotropy, particularly for 
the Heisenberg model. 41-47 While ED is always limited 
to small systems and short times (out to t ~ 6 in units of 
inverse exchange constant for L = 16 sites), 45 QMC 46 is 
plagued by the analytic continuation and cannot resolve 
singularities associated with the long-time behavior. Al- 
though the more recently developed density matrix renor- 
malization group (DMRG) method applied to the trans- 
fer matrix 48 works directly in the thermodynamic limit, 
it is also restricted to intermediate times and has not 
detected a diffusive contribution at low temperatures. 47 
Nonetheless, these works have concluded in favor of the 
existence of diffusion for the Heisenberg model at high 
temperatures. 

Although diffusion was originally proposed to de- 
scribe spin dynamics at high temperatures, the paradigm 
has been used to interpret nuclear magnetic resonance 
(NMR) experiments that measure the spin-lattice relax- 
ation rate 1/T X of spin chains at low temperatures. 49-52 
Strictly speaking, linear response theory expresses 1/Ti 
in terms of the Fourier transform of the transverse spin- 
spin correlation function (Sj~(t)Sj,(t)). The small mag- 
netic field applied in NMR experiments breaks the rota- 
tional spin invariance of Heisenberg chains from SU(2) 
down to U(l) and the total S ± are not conserved. How- 
ever, if the experiments are in the regime of temper- 
atures small compared to the exchange constant, but 
large compared to the nuclear or electronic Larmor fre- 
quencies then the transverse correlation function can be 
traded for the longitudinal one calculated for the elec- 
tronic Larmor frequency. We will discuss this point in 
more detail in Section II G. Up to g-dependent form fac- 
tors that stem from the spatial dependence of the hy- 
perfine couplings, 1/Ti is then proportional to the dy- 
namical autocorrelation G(r — 0,cj), with oj oc h equal 
to the Larmor frequency of an electron in a magnetic 
field h. If diffusion is present, with G(r = 0, t) ~ tr d / 2 
given as in the phenomenological theory, 1/Ti behaves 
as 1/Ti oc G(r = 0,oj) <~ cj d / 2-1 in d dimensions. In 
particular, in the one-dimensional case 1/Ti diverges at 
low frequencies as 1/Ti ~ l/v^ ~ l/y/h. This type 
of behavior has been observed for gapped S = 1 spin 
chains 53 and gapless spin chains with large half-integer 
S. 49 More surprisingly, spin diffusion has also been ob- 
served in S = 1/2 chain compounds by NMR 51,52 and 
by muon spin resonance. 5 It is important to note that 
in the NMR experiment of Ref. [51], the g-dependence of 
the form factor suppresses the contribution from q ~ 7r 
modes in the autocorrelation function, so that the 1/Ti 
signal is completely dominated by q ~ modes. 

The observation of spin diffusion in S = 1/2 Heisen- 
berg chains is puzzling from the point of view of the 
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integrability-transport conjecture. Since the experimen- 
tal diffusion constant was found to be fairly large, 51 it 
becomes important to determine whether the diffusion 
constant is mainly determined by integrability-breaking 
interactions present in the real system or by umklapp 
processes already contained in the integrable Heisenberg 
model. 

Recently, we have shown using a field theory approach 
that the long-time behavior of the autocorrelation func- 
tion and the transport properties are directly related and 
can be obtained from the same retarded Green's func- 
tion at low temperatures. 19 The analytical results were 
supported by DMRG calculations for the time-dependent 
current-current correlation function as well as by a com- 
parison with the NMR experiment on Sr 2 Cu0 3 . 51 We 
also argued that ballistic transport can be reconciled with 
diffusion in the autocorrelation function because ballistic 
channels of propagation can coexist with diffusive ones. 
This can be made precise with the help of the mem- 
ory matrix approach, 10 which allows one to incorporate 
known conservation laws in the low energy effective the- 
ory. However, ballistic and diffusive channels compete 
for spectral weight of the spin-spin correlation function. 
Our field theoretical results for the XXZ chain at h = 
- valid at finite temperatures small compared to the ex- 
change energy - are in very good agreement with the 
diffusive response measured experimentally 51 and with 
time-dependent DMRG results if we assume that the 
Drude weight vanishes completely. Although a small 
Drude weight at finite temperatures cannot be excluded, 
a combination of the numerical data with the memory 
matrix approach implies that it has to be smaller than the 
values obtained in the BA calculation by Kliimper et al. 
[9] and by QMC calculations. 11 - 12 The results in Ref. [19] 
were further supported by a recent QMC study. 30 In the 
latter work the problems arising from analytical continu- 
ation of numerical data were circumvented by comparing 
with the field theory result 19 transformed to imaginary 
times. 

The purpose of this paper is to provide details of 
the calculations for the XX Z chain in Ref. [19]. Fur- 
thermore, we present an extension of these methods to 
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Here Je is the energy current operator as obtained from 
the continuity equation for the Hamiltonian density for 
h = 0. 7,18 It can be verified that Je is conserved in the 
strong sense that [Je,H] = 0. The thermal conductiv- 
ity therefore only has a Drude part which can be calcu- 
lated exactly by BA. 18 According to Mazur's inequality, 
Eq. (1.3), the overlap of J with Je provides a lower 



charge transport in the attractive Hubbard model as 
well as a discussion of the transport properties of the 
Haldane-Shastry chain. Our paper is organized as fol- 
lows: In Sec. II we study the spin current in the spin- 
1/2 Heisenberg chain. We discuss the relation between 
the current-current and the spin-spin correlation func- 
tion at low temperatures, explain in detail how our re- 
sults relate to previous BA and QMC calculations, and 
discuss consequences for electron spin resonance and the 
finite-temperature broadening of the dynamic spin struc- 
ture factor. In Sec. Ill, we discuss spin transport in the 
Haldane-Shastry model and point out that the current 
operator is a nonlocal conserved quantity in the ther- 
modynamic limit. In Sec. IV we show that many of the 
results we obtained for the spin current in the Heisenberg 
model also directly apply to the charge current in the at- 
tractive Hubbard model. Finally, we give a summary and 
some conclusions in Sec. V. 



II. THE SPIN CURRENT IN THE SPIN-1/2 
XXZ MODEL 

The XXZ model (1.4) is exactly solvable by Bethe 
ansatz (BA) 55 and for h = the excitation spectrum is 
gapless for |A| < 1 and gapped for |A| > 1. 

The spin-current density operator is defined from the 
continuity equation for the density of the globally con- 
served spin component 

a t Si = -i[St,H\ = -{j l -j l - 1 ), (2.1) 

which for the XXZ model yields 

Ji = - l -^(StS- +1 -S+ +1 S-). (2.2) 

For h ^ 0, the summed current operator J — ji has 
a finite overlap with the local conserved quantities of the 
XXZ model. The simplest nontrivial conserved quantity 
is 



- SUS?S? +1 ) + A(SUSfSf +1 Sf^SfSfa)] . (2.3) 
I 

bound for the Drude weight of the spin conductivity 

The advantage of this formula is that - contrary to 
Eq. (1.1) - it does not require the calculation of dynami- 
cal correlation functions and is thus much more accessible 
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by standard techniques. The evaluation of (2.4) becomes 
particularly simple in the limit T — > oo leading to 7 
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(T»J), (2.5) 



where m = (Sf) is the magnetization. At low temper- 
atures, on the other hand, standard bosonization tech- 
niques can be applied. Furthermore, BA can be used to 
evaluate (2.4) for all temperatures as will be shown in 
Sec. II A. 

All other local conserved quantities can be obtained ei- 
ther recursively by applying the so-called boost operator 7 
or by taking higher order derivatives of the quantum 
transfer matrix of the XXZ Hamiltonian with respect 
to the spectral parameter. The local operators obtained 
this way act on more and more adjacent sites but they 
are all even under particle-hole transformations. As a 
result, the Mazur bound for the Drude weight, Eq. (1.3), 
vanishes for h — 0. 



Low energy effective model and the Mazur 
bound 



Bosonization of the XXZ model, Eq. (1.4), in the 
gapless regime at zero field leads to the effective 
Hamiltonian 55-57 



H 

H u 
Hbc 



H o + H u + H bc , 

V -Jdx [n 2 + (d x <f>) 2 } , 

A J dxcos(V8irK(f)), 

J dx{dM 2 {d x <t> L ) 2 
J dx [(d x <f> R ) 4 + (d x <t> L ) 4 



(2.6) 



= — 2irv\ + 
- 2irv\- 



Here, H is the standard Luttinger model and H u and 
i?bc are the leading irrelevant perturbations due to Umk- 
lapp scattering and band curvature, respectively. The 
bosonic field <j> — 4>r+4>l and its conjugate momentum II 
obey the canonical commutation relation \(fi(x),Il(x')} = 
i5(x — x'). The long-wavelength (q <~ 0) fluctuation 
part of the spin density is related to the bosonic field by 
SI ~ y/K/2nd x 4>. The spin velocity and the Luttinger 
parameter K are known exactly from Bethe ansatz 
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In this notation, K — 2 at the free fermion point (A = 0) 
and K — 1 at the isotropic point (A = 1). The ampli- 
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In the gapless phase, the Mazur bound can be calcu- 
lated in the low-temperature regime using the field theory 
representations of J and Je- In the continuum limit, the 
continuity equation becomes 

d t S z (x) + d xJ (x) = 0. 

Using the bosonized form for the spin density, we obtain 
for the effective model (2.6) 




dx dt4> 



(2.9) 
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In the following we neglect the corrections to the current 
operator due to band curvature terms and calculate the 
Mazur bound Z?M a zur f° r the Luttinger model H using 
the current operator 



and the energy current operator 

J E w -v 2 J dxnd x 4>. 



(2.10) 



(2.11) 



At zero field, the overlap vanishes because J and Je 
have opposite signatures under the particle-hole trans- 
formation <j) -4- — <^,n — > — LT. A small magnetic field 
term in Eq. (1.4), on the other hand, can be absorbed by 
shifting the bosonic field (here we set \ib — 1) 



(2.12) 
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In this case, the conserved quantity becomes 10 
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= Je + hj. 



(2.13) 

We calculate the equal time correlations within the Lut- 
tinger model and find 
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FIG. 1: Mazur bound -Diviazur as defined in (2.4) for A = 
cos(-7r/4) and magnetic fields h as indicated on the plot. The 
symbols correspond to the exact calculation using the Bethe 
ansatz, the lines represent the field theoretical formula (2.14). 

We note that in the limit T/h — >• the Mazur bound 
obtained from the overlap with Je saturates the exact 
zero temperature Drude weight D(T = 0) = vK/Aw. 32 

One can also use the Bethe ansatz to calculate the 
Mazur bound in Eq. (2.4) exactly. To do so we com- 
puted the equal time correlations using a numerical solu- 
tion of the nonlinear integral equations obtained within 
the Bethe ansatz formalism of Refs. [58,59]. In Fig. 1. 
the numerical Bethe ansatz solution for small values of 
magnetization m is compared to the field theory formula 
(2.14). Remarkably, the free boson result in Eq. (2.14) 
fits well the behavior of the exact Mazur bound out to 
temperatures of order J. 



from decaying. The parameter-free result derived here 
is expected to be correct if the Drude weight is indeed 
zero and allows to check the validity of the assumption 
by comparing with experimental and numerical results. 
Importantly, we can also systematically study how this 
result is modified if such conservation laws do exist after 
all. This case will be considered in section II C. 

We are interested in the retarded spin-spin correlation 
function x ret (<?,u;), which can be obtained from the Mat- 
subara correlation function 

X(q,iu> n ) = ) / dre^ T {St{r)Sm) 

i,v 

(2.15) 

by the analytic continuation iu> n — > lj + i0 + . We will 
show that in the low-temperature limit this correlation 
function determines both the decay of the current-current 
correlation function as well as the spin-lattice relaxation 
rate \jT\. In the low-energy limit, we follow Ref. [60] and 
relate the long-wavelength part of the retarded spin-spin 
correlation function to the boson propagator 

Xretfew) , . rct . . V 

(2.16) 

We calculate the self-energy IF ct (<7, ui) = IF u ot (<7,u;) + 
IK^(q, uS) by perturbation theory to second order in H u 
and first order in H^ c . We first focus on the half-filling 
case (h = 0); the finite field case is discussed at the end 
of this sub-section. 

The contribution from Umklapp scattering reads 

Il™\q,u>) = AkKv\ 2 [F^{ qi u) - F rct (0,0)] , (2.17) 

where 61 



B. Retarded spin-spin correlation function 

We will now proceed with the field theory calculations 
in the following way. We first assume that the Drude 
weight is not affected by unknown nonlocal conserved 
quantities. In this case we are left with the following 
picture based on the Mazur bound and including all con- 
served quantities related to integrability: For finite mag- 
netic field the Drude weight is a continuous function of 
temperature. At zero field, however, the Drude weight 
is only finite at T = but drops abruptly to zero for 
arbitrarily small temperatures. Within the effective field 
theory such a possible broadening of the delta-function 
peak at finite temperatures has to be related to inelastic 
scattering between the bosons which can relax the mo- 
mentum. Such a process is described by the Umklapp 
term in Eq. (2.6) which we have ignored so far. We will 
now include this term as well as the band curvature terms 
in a lowest order perturbative calculation. We want to 
stress that such an approach does not know about con- 
servation laws which could protect a part of the current 



F™\q,w) 



x I 



V 



nT 
v 
vq 



4K 



with 



I(z) = 



du 



2T 



■>2K-1 



sm(2nK) (2.18) 
~2T~ 



o smh 2K (iru) 



-BIK--,1-2K 



where B(x, y) = T(x)T(y)/T(x + y) is the beta function. 
For K > 1/2, we need a cutoff in the integral I(z) in 
Eq. (2.19). However, the imaginary part of I(z) does not 
depend on the cutoff scheme used. 61 The expansion of 
Eq. (2.19) for |w±ug| <C T yields both a real and an imag- 
inary part for IIJ 1 ct (g,w). The calculation of 11^°* (g, uj) is 
also standard. In contrast to IF u ot (<7,u;), the result for 
IL™*(q,oj) is purely real, as band curvature terms do not 
contribute to the decay rate. The end result is 



n ret (<2, oj) « -2i 7 cj - buj 2 + cv 2 q 2 . 



(2.20) 
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In the anisotropic case, — 1 < A < 1, the parameters are 
given by 

n4K-3 



(2.21) 



2 7 = Y{I 
b = (Y 2 -Y 3 )T 4K - 4 + nT^ 

b 2 bi 

c = -(Y 2 + Y 3 )T 4K - 4 -Y 4 T 2 . 

Here b\ and c\ (b 2 and c 2 ) are the parts stemming from 
the band curvature (umklapp) terms, respectively. In 
Eq. (2.21) we have used the following functions 



Y 1 

Y 2 
Y 3 = A 

n = I 



k B(K,l-2K) . T ^ 



B(K,1-2K) 



1 



Tr2 4K + 4 



cot 2 (7rif)r(l/2- K)T(K), 



(2.22) 



r(A++6A_), 



A = AttKX 2 sm{2wK) 



2ir 



AK-2 



T(l/2-K)T(K), 



with ty(x) being the Digamma function. 

At the isotropic point, A = 1, Umklapp scattering 
becomes marginal and can be taken into account by re- 
placing the Luttinger parameter by a running coupling 
constant, K — > 1 + g(T)/2. In this case we find 

4 32 ^ 3 

" 2 3 5 3 VS. 



27 
b 



V3 



T 



(2.23) 



c = 



32 



Following Lukyanov, 57 the running coupling constant 
g(T) is determined by the equation 



In. g 

2 



In 



Tr e 1 /^ 
2 T 



(2.24) 



where 7 is the Euler constant. We remark that a simi- 
lar calculation was attempted in Ref. [13], but there the 
imaginary part of the self-energy was neglected. 

This calculation can be extended to finite magnetic 
field. Shifting the field (f> as in Eq. (2.12), the Umklapp 
term in Eq. (2.6) becomes H u = \J dx cos[V '87T K<p + 
(2Kh/v)x\. As long as h -C T, it is reasonable to keep 
this oscillating term in the effective Hamiltonian. How- 
ever, in a renormalization group treatment, after we lower 
our momentum cut-off below h, it should be dropped. 
Our formula for the self-energy, in Eq. (2.17) is thus mod- 
ified to 



n r u ct (9^) 



= 2nKv\ 2 [F rct (q- 
+ F YCt (q-2Kh/v, 
F ^{2Kh/v,Q)- 



■2Kh/v,w) 
0) (2.25) 
F rct (-2Kh/v,0)] , 



where F rct (q,uj) is given by Eqs. (2.18) and (2.19) as be- 
fore. As a consequence, the relaxation rate for h/T <C 1 
will now be given by 7 - T 4K - 3 [l + 0[(h/T) 2 } + ■ ■ ■} (up 
to logarithmic corrections in the isotropic case). In the 
next section we show that the retarded current-current 
correlation function can be obtained at low energies using 
the calculated self-energy. We will also discuss what the 
shortcomings of the self-energy approach are and show 
that these shortcomings can be addressed by taking con- 
servation laws into account explicitly. 



C. Decay of the current-current correlation 
function 

The time-dependent current-current correlation func- 
tion can be written as 



C(t) 



1 



(J(t)J) 



(2.26) 



f°° du e- lut T l/r ^ , . 



where (J\ J) Te t{w) is the retarded current-current corre- 
lation function. The latter appears in the Kubo formula 
for the optical conductivity (1.1) with (Skin) / L = Kv/2ir 
and the current operator given by J = — yj K /2irdt(t>- 
One can easily show that 

leading to 



v + u 2 {<t>4>) ret {q,u) (2.27) 



Kv K 

{J;JY et {q,u) = -^ + ^cJ 2 {H>) ret {q,u) . (2.28) 
The Kubo formula (1.1) can therefore also be written as 



a(q,uj) 



2tt 



(2.29) 



allowing us to use the results for the boson-boson Green's 
function from the previous section. At zero temperature 
the irrelevant operators in (2.6) can be ignored and the 
Drude weight of the XX Z model can be obtained using 
the free boson propagator 



v 2 q 2 



(2.30) 



implying 

D(T = 0)S{u) 



— lim lim a'(q, tj) 



(2.31) 



Kv i 
-^Re 



Kv 

47T 



5( U ) 



in agreement with Bethe ansatz. 32 

If we now turn to finite temperatures we can use the 
result from the self-energy approach in Eqs. (2.16, 2.20) 
and relation (2.29), leading to the optical conductivity 



a(q,cj) = 



Kv 



2tt (1 + b)uj 2 - (1 + c)v 2 q 2 + 2i 7 w 



(2.32) 
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with the real part being given by 

KvixJ 2-fLJ 



2tt [(l + b)0J 2 -{l + c)v 2 q 2 ] 2 + (2 1 0j) 2 ' 

(2.33) 

For g = 0we find, in particular, a Lorentzian 



vK 



2 7 



2tt [(1 + 6)w] 2 + (2 7 ) 



2 ' 



(2.34) 



As expected, the self-energy approach predicts zero 
Drude weight whenever 7 (T) is nonzero. This result is 
inconsistent with Mazur's inequality if there exist conser- 
vation laws which protect the Drude weight. We know 
that this is the case for any arbitrarily small magnetic 
field h, while our calculations in sub-section IIB gave a 
non-zero 7 at non-zero field. Whether or not such a con- 
servation law exists also for h = is an open question. 
We will try to tackle this problem by first studying how 
the self-energy result is modified by a conservation law, 
followed by a comparison with numerical and experimen- 
tal results. 

It is possible to accommodate the existence of non- 
trivial conservation laws by resorting to the memory ma- 
trix formalism of Ref. [10]. This approach starts from 
the Kubo formula for a conductivity matrix a which in- 
cludes not only the current operator J = J7i, but also 
"slow modes" J n (n > 2) which have a finite overlap with 
J. The idea is that if (J n {t)Jn{fy) is a slowly decaying 
function of time, the projection of J into J n governs 
the long-time behavior of the current-current correlation 
function and consequently dominates the low-frequency 
transport. The overlap between J and the slow modes is 
captured by the off-diagonal elements of the conductiv- 
ity matrix. In practice, only a small number of conserved 
quantities is included in the set of slow modes, but the 
approach can be systematically improved since adding 
more conserved quantities increases the lower bound for 
the conductivity. 62 It is convenient to invert the Kubo 
formula for the conductivity matrix using the projection 
operator method. 62 We introduce the scalar product be- 
tween two operators A and B in the Liouville space 



(A\B) = y dT{A^e HT Be 
L Jo 



Ht 



(2.35) 



Here, C is the Liouville superoperator defined by CJ n = 
[H,J n \. For simplicity, we assumed that all the slow 
modes have the same signature under time-reversal sym- 
metry. The conductivity matrix can be written as 



<T n m(u) = i (J n \{u - £) 1 \J m ) ■ 



(2.36) 



The conductivity in Eq. (1.1) is the an component of a. 
The susceptibility matrix \ is defined as 



Xnm ^ {.\J ri\<Jrn) • 



We denote by 



V — \ T ^""^ Xnm\3n){3rr, 



(2.37) 



(2.38) 



the projector out of the subspace of slow modes. Using 
identities for the projection operator, Eq. (2.36) can be 
brought into the form 10 



^mH=«{[w-Mx *] 1 X}nm, 

where M is the memory matrix given by 



(2.39) 



M nm =T-\j n \LV 



1 



uj-rcr 



PC\J m ). (2.40) 



It can be shown that if there is an exact conservation law 
(local or nonlocal) involving one of the slow modes, the 
memory matrix has a vanishing eigenvalue, which then 
implies a finite Drude weight. 

In the following we apply the memory matrix formal- 
ism to calculate the conductivity for the low-energy effec- 
tive model (2.6) at ft- = 0, allowing for the existence of a 
single conserved quantity Q, [Q,H] = 0. The conductiv- 
ity matrix is then two-dimensional. We choose J\ = J 
and J 2 = Q± = Q- J{J\Q)(J\JY X so that X is diago- 
nal. At low temperatures T <C J, we can use the current 
operator in Eq. (2.10); to first order in A±, we obtain 



X11 



(J 2 ) vK 



LT 



2tt 



(2.41) 



where b\ is defined in Eq. (2.21) and we have used the 
fact that due to the vanishing of the superfluid den- 
sity we have (E^^/L w (J 2 )/LT. Likewise, X22 = 
((Q 2 ) - (QJ) 2 /(J 2 ))/(LT) can be calculated within the 
low energy effective model once a conserved quantity Q 
has been identified. The remaining approximation is in 
the calculation of the memory matrix to second order in 
Umklapp. This is analogous to the calculation of the self- 
energy IT u ct in Eq. (2.17). Using the conservation law, we 
can write 



M w M11 (a;) 



where r = (QJ)/(J 2 ) and 



1 - 
-r r 



(2.42) 



MnH^^^U.-^), (2.43) 

Z7T Ul 111 

with 6 2 as given in Eq. (2.21). Thus, from equation 
(2.39), we find 



o-(w) = X11 
where 

V = 



y i 1 

- + 



l + yuj l + yco-{l + y)xuM n {oj)_ 

(2.44) 



r 2 Xn 



(JQY 



X22 {J 2 ){Q 2 )-{JQY 



(2.45) 



For finite magnetic field, a conserved quantity is given by 
Q = J E + hJ (see Eq. (2.13)) and thus y - (h/T) 2 in 
this case. Equating Eq. (2.32) for q — and (2.44), we 
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find that the memory matrix approach is equivalent to 
adopting the self-energy 



n( w ) 



l-yXn-WnM/w 



- - yXn^iiM/^ 

(2.46) 

As expected, the memory matrix result reduces to the 
self-energy result for y — > 0. The difference between the 
self-energy and the memory matrix result is of higher 
order in the Umklapp interaction. Therefore, the con- 
servation law is not manifested in the lowest-order cal- 
culation of the self-energy. Although Eq. (2.46) is not 
correct beyond (9(A 2 ,A±), it suggests that the memory 
matrix approach corresponds to a partial resummation 
of an infinite family of Feynman diagrams which changes 
the behavior of n(w) in the limit lu — > from II — > —2ijuj 
to II -)• -jr 1 ^ 2 . 

It follows from Eq. (1.1) and Eq. (2.44) that 



Im{J;J) Ict (oj) = -wRea(w) 



vK 



+ 



7[j/(l - frl) 

l + y 

27 



S(oj) (2.47) 



(! + &! + 6' 2 ) 2 ^ 2 + (2 7 ') 2 



where b 2 = (1 + y)&2 and 7' = (1 + y)-f. The first term 
on the right hand side of Eq. (2.47) can be associated 
with the ballistic channel and the second one with the 
diffusive channel. The calculation is valid also in the fi- 
nite field case if h/T <C 1 with y <~ (h/T) 2 as already 
discussed below Eq. (2.44). This means that we obtain 
a correction of the relaxation rate in the memory matrix 
formalism which is exactly of the same order as obtained 
previously from the self-energy approach in this limit (see 
Eq. (2.25)). However, now we see that the weight is trans- 
ferred accordingly from the diffusive into a ballistic chan- 
nel, a fact which is missed in the self-energy approach. 

Substituting Eq. (2.47) into Eq. (2.26) we can now 
evaluate the integral. From the second term in (2.47) 
we obtain an integrand which has poles at frequencies 
uj = ±2i7'/(l + h + b 2 ) and w = 2mTn with n e Z. 
Since 27' <C 2txT the contributions of the latter poles 
can be ignored at times t >• (2nT)^ 1 leading to 



C(t) 



Kv 



2tt(1 + y) 



yT(i - h) 



2i 7 " 



-2~/"t 



1 + h + b' 2 1 - e 2 n"/T 
(2.48) 

with 7" = 7 ; /(l + b\ + b' 2 ). We can further expand the 
denominator in powers of 7" /T <C 1. The first order 
contribution is real and given by 



C(t) 



vKT 
2tt(1 + y) 



y(l-6i) + 



-2Vt 



1 + 61 + b' 2 



(2.49) 



The imaginary part of the correlation function is obtained 
in second order in the expansion and is thus suppressed 
by an additional power of 7 '/T . From Eq. (2.49) we see 



that in the limit t — > 00 the current-current correlation 
function approaches the value 

= vKTyjl-h) ^ (J 2 )y = (JQ) 2 
t^- [) 2n(l + y) ~L(l + ») L(Q 2 ) ' 

(2.50) 

consistent with the Mazur bound for the Drude weight. 
For intermediate times (2ttT)^ 1 <C t <C I/7', we obtain 
the linear decay 



C(t) 



vKT 
2tt(1 + b) 



(1 - 2t*), 



(2.51) 



independent of y if 61, 6 2 ^ 1- Therefore a small 
Drude weight cannot be detected in this intermediate 
time range. 



D. Comparison with Bethe ansatz and numerical 
data for the current-current correlation function 

In the previous section we have derived a result for 
the optical conductivity, Eq. (2.47), and for the real part 
of the current-current correlation function, Eq. (2.49), 
at low temperatures by a memory-matrix formalism. In 
this approach we have explicitly taken into account the 
possibility of a nonlocal conservation law - leading to a 
nonzero Drude weight at finite temperatures. By com- 
paring our results with the Bethe ansatz calculations 8 ' 9 
and numerical calculations of C(t) we will show that a 
diffusive channel for transport does indeed exist. Fur- 
thermore, we will try to obtain a rough bound for how 
large y and therefore the Drude weight D(T) can possibly 
be. 

A finite Drude weight at finite temperatures has 
been obtained in two independent Bethe ansatz 
calculations, 8 ' 9 however, the obtained temperature de- 
pendence is rather different. In both works the finite 
temperature Kohn formula, 6 which relates the Drude 
weight and the curvature of energy levels with respect 
to a twist in the boundary conditions, is used. While 
Zotos [8] uses a TBA approach based on magnons and 
their bound states, Klumper et al. [9] use an approach 
based on a spinon and anti-spinon particle basis. Both 
approaches have been shown to violate exact relations 
at high temperatures and are therefore not exact solu- 
tions of the problem. The reason is that within the Bethe 
ansatz both approaches use assumptions which have been 
shown to work in the thermodynamic limit for the par- 
tition function. This, however, does not seem to be the 
case for the curvature of energy levels relevant for the 
Drude weight. Nevertheless, this does not exclude that 
these results become asymptotically exact at low temper- 
atures and arguments for such a scenario have been given 
in Ref. [9]. 

To investigate this possibility we start by comparing 
in Fig. 2 the Bethe ansatz results from Ref. [9] at low 
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0.16r 



0.159 




0.156 



0.02 0.04 
T/J 



0.02 0.04 0.06 
T/J 



FIG. 2: The Bethe ansatz results from Ref. [9] (dots) com- 
pared to the field theory formula (2.52) (lines) obtained by 
setting 7 = by hand in Eq. (2.32). The A values are indi- 
cated on the plot. 



temperatures with the Drude weight 



D(T) = 



Kv 



4tt(1 + b) 



(2.52) 



obtained by setting 7 = in Eq. (2.32). We see from 
Eq. (2.49) that this corresponds to the case y — > 00. i.e.. 
in this case there is only a ballistic channel. Doing so we 
obtain excellent agreement. From this we draw two con- 
clusions: First, this BA calculation predicts that even at 
finite temperatures the transport is purely ballistic. Sec- 
ond, the temperature dependent parameter b which we 
obtained from field theory in first order in band curva- 
ture and second order in Umklapp scattering is consistent 
with this BA approach. We note that the case A = 1/2 
(K = 3/2) is particularly interesting because here the 
contributions b\ from band curvature and 6 2 from Umk- 
lapp scattering in Eq. (2.21) both yield a T 2 contribu- 
tion with diverging prefactors. These divergencies cancel 
leading to a Drude weight 



D(T) 



9%/3 



1 



32tt l + bT 2 



(2.53) 



with b = [142+24 7 +48 ln2+60 In 3-241nT-21C(3)]/243 
where 7 is the Euler constant and £ the Riemann zeta 
function. As is also the case for other thermodynamic 
quantities 63 we see that the coinciding scaling dimen- 
sions of the terms stemming from Umklapp scattering 
and band curvature lead to a term ~ T 2 In T at this 
special point. For A < 1/2 (K > 3/2) band curva- 
ture gives the dominant temperature dependence while 
Umklapp scattering dominates for A > 1/2 (A" < 3/2). 
For K > 5/4 (A > 0.81) a T SK ~ S term, which arises 
in 4th order perturbation theory in Umklapp scattering 
and which is not included in our calculations, becomes 
more important than the T 2 term from band curvature. 



Therefore the agreement for A = 0.9 in Fig. 2 is not 
quite as good as for the other values. For A = 1, Umk- 
lapp scattering becomes marginal leading to a logarith- 
mic temperature dependence of the Drude weight (2.52) 
with b as given in (2.23). 

To see whether or not such a large Drude weight as pre- 
dicted by Khimper et al. is possible and to discuss the 
second BA approach by Zotos we now turn to a numeri- 
cal calculation of C(t) . A dynamical correlation function 
at finite temperatures can be obtained by using a den- 
sity matrix renormalization group algorithm applied to 
transfer matrices (TMRG). 19 ' 47 * 48 This algorithm uses a 
Trotter-Suzuki decomposition to map the ID quantum 
model onto a 2D classical model. For the classical model 
a so-called quantum transfer matrix can be defined which 
evolves along the spatial direction and allows one to per- 
form the thermodynamic limit exactly. In order to cal- 
culate dynamical quantities, a complex quantum transfer 
matrix is considered with one part representing the ther- 
mal density matrix and the other part the unitary time 
evolution operator. By extending the transfer matrix one 
can either lower the temperature (imaginary time) or in- 
crease the real time interval for the correlation function. 
The calculation of the current-current correlation func- 
tion is particularly complicated because it involves the 
summation of all the local time-dependent correlations 



5>l(*)io(0)) (2.54) 



with the current density ji as given in Eq. (2.2). In Fig. 3 
the local correlations (ji(t)jo(0)) are exemplarily shown 
for the case A = 0.6 and T/J — 0.2. In order to obtain 




FIG. 3: (ji(t)jo(0)) for A = 0.6 at T/J = 0.2. 

converged results, the two-point correlations for distances 
up to I 30 have to be summed up. This is not possible 
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using exact diagonalization which is restricted to consid- 
erably smaller system sizes. A good check is obtained by 
considering the free fermion case A = 0. Here (ji(t)jo(0)) 
can be calculated exactly and is non-trivial, however, 
J = ii commutes with the Hamiltonian leading to 
(2.54) being a constant. 

We now discuss the same parameter set A = 0.6 and 
T/J = 0.2 as used above in more detail. In Fig. 4 the real 
and imaginary parts of C(T)/2JT, obtained in a TMRG 
calculation with a 1000 states per block kept, are shown. 
Because the imaginary part is very small, an extrapola- 
tion in the Trotter parameter was necessary restricting 
the calculations to smaller times than for the real part. 
In the limit t — > oo, the real part would directly yield the 




0.155 

H 




-0.015, 



FIG. 4: Real and imaginary parts of the current correlation 
function for A = 0.6 and T/J = 0.2. The dashed, dot-dashed 
and dotted lines correspond to formula (2.48) with y values 
as indicated on the plot. 



Drude weight (possibly zero). As already discussed, the 
BA result by Kliimper et al. corresponds to purely ballis- 
tic transport, y — > oo. From Fig. 4 we see that this is not 
consistent with the numerical data. The BA calculation 
by Zotos, on the other hand, predicts a Drude weight 
~ 0.105 which requires y ~ 2.5. While formula (2.48) 
with y = seems to fit the numerical results best, we 
would need to be able to simulate slightly longer times 
(a factor of 1.5 — 2 should be sufficient) to clearly dis- 
tinguish between y = and y = 2.5. Most importantly, 
however, the numerical data clearly demonstrate that the 
decay rate 7 is nonzero. 

In addition to the zero field case, we have also calcu- 
lated C(t) at relatively large magnetic fields and various 
temperatures. As shown in Fig. 5 we find that in such 
cases C(t) appears to converge to a finite value within 
fairly short times. Furthermore, we find (see Fig. 6) that 
at large A and large magnetic fields the Mazur bound 
(2.5) almost completely exhausts the Drude weight at 
high temperatures. This is consistent with the findings 
in Ref. [7] which were based on exact diagonalization. 



FIG. 5: C(t) for A = 0.4 and a magnetization (m) = 0.3 
at temperatures T/J = 20,1,0.6,0.4,0.2 (in arrow direc- 
tion). At low temperatures the asymptotic value seems to 
be reached almost instantaneously. At high temperatures 
the data seem to be consistent with a simple exponential 
decay to a finite value without oscillations. Using an expo- 
nential fit and extrapolating to infinite temperature we find 
C(T -> 00, t ->• 00) « 0.058. On the other hand, the Mazur 
bound, Eq. (2.5), yields 2TL> Mazur = 0.013. 



0.08 



U 




0.06 



0.04 



0.02. 



FIG. 6: C(t) for A = 4.0 and a magnetization (m) = 0.25 
with T/J — 1, 5, 10, 20 (in arrow direction). In addition, lin- 
ear fits of the data in the regime Jt > 4 are shown (dashed 
lines). Extrapolating these values we find C(T — > 00, t — > 
00) w 0.0385. In this case the Mazur bound, Eq. (2.5), yields 
2T - D Mazur = 0.0366 for T -> 00 and exhausts 95% of the 
Drude weight. 



E. Comparison with Quantum Monte Carlo 

Quantum Monte Carlo (QMC) calculations are per- 
formed in an imaginary time framework. In Refs. [11,12] 
the optical conductivity a(q, iu) n ) at Matsubara frequen- 
cies w„ = 27rTn, n € N has been determined. In order to 
answer the question whether or not a finite Drude weight 
exists at finite temperatures, one has to perform first the 
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FIG. 7: QMC data (symbols) for T/J 

L = 128 for (a) A = 1/2, and (b) A = l. 04 In comparison 
the field theory result in Eq. (2.32) (solid lines) and the result 
obtained by setting 7 = (dashed lines) is shown. For A = 
1/2 both cases are almost indistinguishable. For A = 1 and 
small q the effect of 7 is largest, however, in this case also 
the finite size corrections are large so that a detailed analysis 
seems very difficult. 



limit q — ¥ and then try to extrapolate in the discrete 
Matsubara frequencies u ui n —> 0". Doing so the results 
obtained in Refs. [11,12] have been interpreted as being 
consistent with the Drude weight found in the BA calcu- 
lations by Kliimper et al. [9]. 

By replacing uj — > iu> n , the field theory formula (2.32) 
yields a prediction for a(q,iuj n ) which can be directly 
compared with the QMC results (see Fig. 7). In addi- 
tion we also show the result that is obtained by setting 
7 = in Eq. (2.32), corresponding to the BA solution 
in Ref. [9]. Since the agreement is good for both 7 = 
and 7 as given by Eq. (2.21), we conclude that these 
QMC calculations are not of sufficient accuracy to decide 
whether or not the relaxation rate vanishes in the inte- 
grable model. A general problem in QMC calculations is 
that a relaxation rate much smaller than the separation 
between Matsubara frequencies, 7 <C u) n +i ~ w n = 27rT, 
cannot be resolved. 



Promising seems to be a study of the Heisenberg point, 
where "f(T) ~ T/ln 2 (J/T) is largest. Indeed, evidence 
for a nonzero relaxation rate for this case has been found 
in a very recent QMC study by directly comparing with 
our result (2.32) transformed to Matsubara frequencies. 30 
The magnitude of 7 seemed to be roughly consistent with 
the value in Eq. (2.23). However, as has been already 
shown in Ref. [19] by comparing with numerical results 
for the correlation function C(f), logarithmic correction 
at the isotropic point limit the temperature range where 
the field theory results are applicable. In QMC calcula- 
tions a further problem at A = 1 is the slow logarithmic 
decay of finite size corrections. It would therefore be de- 
sirable to perform a similar study for A ~ 0.6 — 0.8 where 
7 is still fairly large and the field theory seems to work 
well up to temperatures of the order T/J ~ 0.2 as we 
demonstrated in the previous section and in Ref. [19]. 



F. Spin diffusion 

So far we have used our result for the self-energy of 
the bosonic propagator to discuss the transport proper- 
ties of the spin chain. In this and the following sections 
we will discuss diffusive properties characterized by the 
long-time behavior of the spin-spin correlation function 
(Si +x (t)Sf(0)). We concentrate again on the case of zero 
field. For T = 0, it is known 65 that the slowest decaying 
term in the autocorrelation function for < A < 1 is of 
the form 



-iWt 



(T = 0) 



(2.55) 



with W — v and 77 = {K + l)/2. This oscillating term 
is attributed to q = Tt/2 high-energy particle-hole exci- 
tations with a hole near the bottom of the band and a 
particle at the Fermi surface, or a particle at the top of 
the band and a hole at the Fermi surface. At T = 0, 
the low-energy contributions to the spin-spin correlation 
function decay faster than the high-energy contributions. 
In contrast, numerical studies seem to suggest that at 
high temperatures the oscillating terms are still present, 
but the slowest decaying term has pure power-law decay 
with no oscillations. 45 ' 47 This slowly decaying term has 
been interpreted as due to spin diffusion at high temper- 
atures. Here we will show that a diffusive term is already 
present at low temperatures. In the following we use 
the boson propagator in Eq. (2.16), (2.20) to calculate 
the q ~ low-energy contribution to the autocorrelation 
function in the regime T« J. 

We can write the low-energy, long-wavelength contri- 
bution to (Sf +x (t)Sf (0)) as 



G{x,t) 



+°° duj_ 
.~ 2^ 



+00 



dq 
2^ 



ImXretfew) 



1 



e - w /T ' 
(2.56) 

with Xrct(ciS w ) given by Eq. (2.16) and (2.20) . Note 
that the — btu 2 and cv 2 q 2 terms in the self-energy sim- 
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ply rescale the energy and wave- vector by (1 + b) and ^/(l + c)/(l + b)v. Doing the integral over q first, we 
(1 + c) respectively, giving a T-dependent velocity, v = find (see appendix A for details) 



G{x,t) 



K 

2-KV' 



(l + 6)-2(i + c )- 



du) 

47T 



exp[-iwt- ^ 2 + 2tju)^ 2 \x\/v] + 2 ^) _ (2.57) 

1 — e^^i 1 



where t = t(l + b)-^ 2 , x = ^(l + c)- 1 / 2 , 7 = 7(l + &)~ 1 / 2 
and f = T(1 + b) 1 ' 2 . The integral in Eq. (2.57) has both 
pole and branch cut contributions, inside the light cone 
when \t\ > v\x\, but only the pole contribution outside 
the light cone when \t\ < v\x\. While this integral could 
be evaluated more generally numerically, we focus on a 
couple of simple regions where we can obtain analytic 
results. One of these is \x\/v <gc 1/7 <C t where 

G{x,t) = (l + b)-^ 2 (l + c)- 3 / 2 [G (x,t)+G int (x,t)}, 

(2.58) 

with 

G (x,t) = I J= } (2.59) 



K 



nT 



Stt 2 v 2 I sinh[7rf(t + x/v)] j 

is the Luttinger liquid result in terms of the rescaled vari- 
ables and 



Gi nt (x,t) 



KT 

1 12 



-~/x 2 /2v 2 t 



(2.60) 



is the diffusive term. The other simple region is \t\ < 
\x\/v <C I/7 where we just obtain the Luttinger liq- 
uid result, Go- Note that even for the static correla- 
tion function there are for \x\/v > 1/7 in principle 7- 
dependent corrections to the approximately exponential 
decay (see appendix A). However, since I/7 ^ £ where 
£ = v/(2ttT) is the correlation length, these corrections 
occur in a regime where the correlation function is al- 
ready extremely small and where also higher order correc- 
tions to (2.59) have to be taken into account. A schematic 
representation of the behavior expected in the different 
regions is shown in Fig. 8. We expect these results to 
be universally valid at large x and t, ignoring a possible 
Drude weight term, since the Fourier transform is domi- 
nated by the small uj and q region, at large t and x. G m t 
has the classic diffusion form. Note that diffusion oc- 
curs in a more limited domain than originally proposed 
for high-T ferromagnets 38-40 however. We have shown 
it to occur at low T, where 7 <§; T, ignoring a possible 
non-zero Drude weight term in the self-energy, but only 
inside the light cone at \x\/v 1/7 <C \t\ and only for 
the uniform part of the spin correlation function. (There 
is also a power-law oscillating term.) Importantly, the 




Light 
cone 



§fr Luttinger 
liquid 



|x|/v 



1/7 



FIG. 8: The low-energy, long- wavelength contributions to 
{Si +x (t)S z (0)} at finite temperatures in different regions of 
the spacetime diagram. Inside the light cone, Eq. (2.57) has 
both branch cut and pole contributions but only pole contri- 
butions outside the light cone. The Luttinger liquid result is 
valid for \t\ < \x\/v <C I/7 whereas the diffusive term (2.60) 
dominates for \x\/v <g£ I/7 <C t. 



domain where classical diffusive behavior occurs does in- 
clude the self-correlation function. For all < A < 1, at 
finite temperatures and sufficiently long times the auto- 
correlation function becomes dominated by a low-energy 
term 



(Sf(t)Sf(0)) ~ T 



7(r) 
t 



(2V0,t»l/7) (2.61) 



with the universal power-law decay i -1 / 2 expected for 
diffusion in one dimension. 38-40 Note that the more lim- 
ited diffusive behavior we have found is related to the 
Lorentz invariance of the underlying low energy effective 
Lagrangian. The b and c terms break Lorentz invari- 
ance but only produce an unimportantT-dependent shift 
of the velocity. The important breaking of Lorentz in- 
variance, which leads to diffusion, is due to the finite 
temperature. While diffusive behavior occurs only in a 
limited space-time domain, it is sufficient to give diffusive 
behavior in certain NMR experiments, as we show in the 
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next sub-section. 



G. Spin-lattice relaxation rate 

If the q ~ contribution dominates the dynamics, 
diffusive behavior should show up as characteristic fre- 
quency and magnetic field dependence in NMR experi- 
ments. Before deriving a prediction for the spin-lattice 
relaxation rate based on the results from the previous 
section and comparing with experiment, we first want to 
point out a general relation useful to take the finite mag- 
netic field in NMR experiments into account properly. 

The linear response formula for the spin-lattice relax- 
ation rate is 66 



rate 



1 



dq 
2^ 



|A(g)| 2 S+-( g , Wjv )L, 



(2.62) 



where A(q) is the hyperfine coupling form factor, ujn is 
the nuclear magnetic resonance frequency and 

1 r +oc 

S + ~(q,u;)\ h = -J2 / dte^(S+(t)S-(0))\ h 

1,1' 

(2.63) 

is the transverse dynamical spin structure factor. Here 
Sf 1 = Sf ± iSf are the raising and lowering spin oper- 
ators. The expression in Eq. (2.63) is to be calculated 
using Hamiltonian (1.4) in the presence of a magnetic 
field h. We focus on the experimentally relevant Hciscn- 
berg point A = 1. We would like to express 1/Ti in terms 
of the longitudinal structure factor 



S**(q,w) 



1,1' J -°° 



dte lUNt (S?(t)S?,(0)) (2.64) 



at zero field. The latter is more easily calculated in the 
field theory since Sf is related to the local density of 
fermions n;, whereas S t have nonlocal representations 
in terms of Jordan- Wigner fermions. Although the ex- 
change term in the Heisenberg model is isotropic, the 
magnetic field term in Eq. (1.4) breaks rotational sym- 
metry. As a result, S A (q,tj)\h cannot be directly re- 
placed by 2S zz {q 1 uj)\h at finite field. However, we note 
that the longitudinal field has a trivial effect on Sf (t) 



jm s ± e -im 



(2.65) 



where H — H(h — 0) and uj e = fish is the electron mag- 
netic resonance frequency. If we assume in addition that 
T ^> w e , the magnetic field dependence in the thermal 
average can be neglected and we have 

2 r+°o 

S + -(q,u)\ h / dte^-^iSn^mio, 

l,V J -°° 

(2.66) 

where the correlation function is calculated at ft, = 0. 
This leads to the expression for the spin-lattice relaxation 



|r«/ ^\A(q)\ 2 S"(q,w N -w e )\o. (2. 



T 



Using oj e ^> ujn and 



S"(q,-w e ) 



21m Xretfa^e) 

1 _ e ^/T ■■ 



67) 



(2.68) 



we find in the regime w e < T 

1 _2T f ^ ]A{q)] 2 Im Xrct((ZjCJe) . ( 2. 6 9) 
±i u! e J 2ir 

If the integral in Eq. (2.69) is dominated at low temper- 
atures, T -C J, by the q ~ mode of Xret(?>w) then we 
can perform the momentum integral using the retarded 
correlation function in Eq. (2.16). Assuming further- 
more that the momentum dependence can be neglected, 
A(q ~ 0) = A = const, and using the appropriate pa- 
rameters (2.23) for the isotropic case we find 



1 

jyf 

where 



2\A\ 2 X 




, (2.70) 



g g 2 , 5 5 3 , 3V3^ 



64 



2tt 



+ —T. (2.71) 

7T 



In the limit j(T) ^> oj e , we obtain the diffusive behavior 



h{T) 



' T I In 2 {J jT) 



(2.72) 



In the NMR experiment by Thurber et ai, Ref. [51], 
the copper-oxygen spin chain compound Sr 2 Cu03 was 
studied. In this compound an in-chain oxygen site exists 
such that A(q) = Acos(q/2) in (2.69). In NMR measure- 
ments on this oxygen site any contribution to the spin- 
lattice relaxation rate coming from q ~ 7r is therefore 
almost completely suppressed so that the q ~ mode 
dominates. Note that for NMR measurements on the 
copper or the apical oxygen site both low-energy contri- 
butions would be present with the q ~ it being dominant. 
In Ref. [51] the experimental data for the in-chain oxy- 
gen site have been interpreted in terms of a spin-lattice 
relaxation rate \jT{F <~ T/y/uJ^. While the frequency 
dependence does agree with that found in Eq. (2.72) and 
generally expected if diffusion holds, we note that the 
temperature dependence of the effective diffusion con- 
stant D s = v 2 /27 is different. 

To quantitatively compare our prediction with experi- 
ment we note that the form factor A(q) = Acos(q/2) for 
the in-chain oxygen site is given by 



\A\ 2 = 



k B (2C b ) 2 + (2C C Y 
2h tt 3 /c|J 2 



(2.73) 
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where ks is the Boltzmann constant, C b,c are the di- 
mensionless components of the hyperfine coupling tensor, 
giN^- = 4.74 x 10 -9 eV and J is the exchange coupling 
measured in Kelvin. The components of the hyperfine 
coupling tensor can be obtained by performing a K — \ 
analysis, i.e., a comparison between measurements of the 
Knight shift K and the magnetic susceptibility \. Such 
an analysis has been performed in Ref. [51] leading to 
2C b = 95 ±10 and 2C C = 44 ±10 kOe/> B . Furthermore, 
the exchange coupling J has to be determined. We note 
that a parameter-free field theory formula for the static 
magnetic susceptibility follows from 



-XretfeO) 



1 



K 1 

2ttv 1 + c 

32 



(2.74) 



where we have used Eqs. (2.16, 2.20). In the second line, 
we have specialized for the isotropic point, with v = tt/2, 
K —> 1 + g/2 + <7 2 /4 + g s /8, and the parameter c as given 
in (2.23). Eq. (2.74) is in agreement with the result in 
Refs. [57,67]. Using this formula and fitting the suscep- 
tibility data in Ref. [68] we obtain J = 2000 ± 200 K. 
In Fig. 9 the shaded area denotes the region covered by 
the curves obtained by varying the parameters J and C b ' c 
independently within the given range. However, if we be- 
lieve that the field theory describes the zero temperature 
limit correctly then the possible variation of the theoreti- 
cal results is strongly overestimated. From Eq. (2.70) we 
see that 1/{T\T) -> 2\A\ 2 for T — > and is therefore not 
affected by the question how large 7 is. For a given J 
we can therefore fix (2C b ) 2 + (2C C ) 2 in (2.73) by fitting 
the experimental data at low temperatures. The remain- 
ing variation is then very small (see the two examples 
in Fig. 9). For 7 = corresponding to purely ballistic 
transport, as suggested by the Bethe ansatz calculation 
of Ref. [9], 1/(T\T) would be approximately constant. 19 
Further support that the good agreement seen in Fig. 9 
is not accidental is obtained by an analysis of the NMR 
data for the apical oxygen site which is also presented 
in Ref. [51]. In this case the hyperfine coupling ten- 
sor A(q) is momentum independent and the contribution 
from low-energy excitations with q ~ ir dominates. In 
this case field theory yields 69 



1 



(D b y + (D c f { giN ny 



(2tt) 3 / 2 



hknJ 



(2.75) 



with D b = 23 ± 10 and D b = 14 ± 10 kOe//z B . 51 Using 
again J = 2000 K the experimental data in Ref. [51] are 
very well described by this formula. Therefore the experi- 
mental data for the magnetic susceptibility and the NMR 
relaxation rates, testing the low-energy contributions of 
Xret(<7>' t ')> both at q ~ and q ~ w, are consistently de- 
scribed using J ~ 2000 K and a hyperfine coupling tensor 
as determined experimentally by a K — \ analysis. 
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FIG. 9: Experimental data for the spin-lattice relaxation rate 
of the spin chain compound Sr2CuC>3 at h = 9 T taken from 
Ref. [51] (dots) compared to our theory. The solid lines rep- 
resent the two extreme limits when varying J, C b ' c indepen- 
dently within the given error bars. The other two curves 
correspond to fixing the hyperfine constants for a given J by 
the experimental data at low temperatures. 



H. Consequences for electron spin resonance 

Here we want to point out a connection between our 
study of the conductivity of the 5=1/2 Heisenberg chain 
and a previous theory of electron spin resonance (ESR) 
for such systems. 60 By also using a self-energy approach 
for the boson propagator it was shown in this work that 
the ESR lineshape is Lorentzian with a width given by 
the imaginary part of the retarded self-energy. In partic- 
ular, the case of an exchange anisotropy perpendicular to 
the applied magnetic field was considered. In this case 
the ESR linewidth for an applied magnetic field h is given 

by 



V = -lmUT(h,h)/(2h) 



(2.76) 



where n* et (ft,, h) is the retarded self-energy due to Umk- 
lapp scattering, Eq. (2.17), for the isotropic case. From 
(2.20) and (2.23) we therefore immediately obtain 



T] EE 7 



;9 2 T. 



(2.77) 



This is consistent with Eq. (6.27) in Ref. [60] using the 
standard replacement A — > g/A in reverse and in addition 
A — > vX with v = tt/2 because in [60] the spin velocity was 
set to 1. The RG flow for the running coupling constant 
g is now cut off by the larger of the temperature T or 
the applied magnetic field h. If these two scales are suffi- 
ciently different, we can use Eq. (2.24) for g with T being 
replaced by max(T, h). In this case (2.77) is a parameter- 
free prediction for the ESR linewidth of a S — 1/2 Heisen- 
berg chain with a small exchange anisotropy perpendic- 
ular to the applied magnetic field. 

We want to stress that a vanishing relaxation rate 
7 would dramatically affect the ESR linewidth. From 
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(2.17) and (2.18) we obtain for u/T, q/T < 1 



III. HALDANE-SHASTRY MODEL 



(2.78) 



2,2, 9 2 ^ 2 , • fr, , 



This leads to 

r? = -Imn^ ct (/i, h)/(2ft) - 7 + /i 2 /T (2.79) 
and therefore a line width 77 ~ /i 2 /T if 7 = 0. 

I. Consequences for the spin structure factor 

The longitudinal dynamical spin structure factor at fi- 
nite temperatures can be obtained from 



S zz (q,u) = - 



1-e- 



-lm X ret {q,uj) 



(2.80) 



In the low-energy, long-wavelength limit the retarded 
spin-spin correlation function can be again expressed by 
the boson propagator using Eq. (2.16). With the help of 
(2.29) we therefore obtain a direct relation between the 
spin structure factor and the real part of the conductivity 

S"{q,u>) = ^ ^l-prf 'M (<?«!)• (2-81) 

>From (2.33) we see that our theory predicts a 
Lorentzian lineshape at finite temperatures. For T — > 
the parameters 6,0,7 m (2.33) vanish and we find the 
well-known free boson result 2S zz (q,uj) = K\q\5(uj—v\q\). 

Here some comments are in order. In our perturba- 
tive calculation we have only included the band curvature 
terms to first order. While the lowest order contribution 
vanishes at zero temperature, the terms of order 2n show 
divergencies on shell, uj ~ vq. 70 While it is not yet clear 
how to sum up a series of such terms, we know from BA 71 
that they lead to a finite linewidth 5ui(q) ~ q 3 at zero 
temperature. The lineshape at zero temperature is dis- 
tinctly non-Lorentzian with threshold singularities. For 
finite temperature such that 7(T) ^> 5uo(q) we can, how- 
ever, neglect these terms and Eq. (2.81) becomes valid. 
Furthermore, we note that the second order Umklapp 
contribution to the self-energy leads for zero temperature 
to the high-frequency tail of S zz {q 1 uj) given in Eq. (7.24) 
of Ref. [70]. From Eq. (2.81) it is also clear that the regu- 
lar part of the conductivity and a part which could possi- 
bly yield a ballistic channel for dc transport are not inde- 
pendent but rather have to fulfill sum rules, as for exam- 
pie, JduS"(q,w) = 2nJ2 J exp(iqj)(S 1 (0)S a (0)). The 
right hand side of this sum rule can be calculated nu- 
merically with high accuracy since it only involves static 
correlation functions. 



Unlike the XX Z model, generic spin models with ex- 
change interactions beyond nearest neighbor are not in- 
tegrable. Here we want to consider a special spin model 
with long-range interactions that is known to be inte- 
grable, even though it is not solvable by Bethe ansatz. 
The Haldane-Shastry model is given by 



H 



3<l 



(3.1) 



where Jji = J(j - I) = {(L/tt) sin[7r(j - l)/L}}~ 2 is the 
long-range exchange interaction. 72 ' 73 For a chain with 
size L and periodic boundary conditions, Jji is inversely 
proportional to the square of the chord distance be- 
tween sites j and I. We set the energy scale of the ex- 
change interaction such that in the thermodynamic limit 
Jji -+ (j 

The Haldane-Shastry model is completely integrable 
as the transfer matrix satisfies the Yang-Baxter equa- 
tion. The conserved quantities are nonlocal and are not 
obtained by the traditional method of expanding the 
transfer matrix. However, they have been found using 
the "freezing trick" starting from the SU(2) Calogero- 
Sutherland model. 74 The exact spectrum is given by the 
dispersion of free spinons, which contrasts with the in- 
teracting spinons of the Heisenberg model. The dynami- 
cal structure factor has been calculated exactly. 75 It has 
a square-root singularity at the lower threshold of the 
two-spinon continuum, as expected for SU(2) symmetric 
models. Contrary to the Heisenberg model, there are no 
multiplicative logarithmic corrections to the asymptotic 
behavior of correlation functions. This means that in the 
effective field theory for the Haldane-Shastry model the 
coupling constant of the umklapp operator is fine tuned 
to zero. 

Since the umklapp operator is responsible for the lead- 
ing decay process of the spin current in the calculation 
in Sec. II C, we may expect that spin transport in the 
Haldane-Shastry model is purely ballistic. This seems 
reasonable given the picture of an ideal spinon gas, but 
requires that the coupling constants of all higher order 
Umklapp processes also vanish exactly. 

Although the total spin is conserved, the spin density 
does not satisfy a local continuity equation due to the 
long-range nature of the exchange interactions. Nonethe- 
less, the spin current operator can be defined using the 
equivalence to a model of spinless fermions coupled to 
an electromagnetic field. In terms of Jordan- Wigner 
fermions, the Haldane-Shastry model reads 



H 



E4« 



3<l 



n j 



+ 



(-!)< 



ni 



h 



,) 



(3.2) 



We then consider a magnetic flux $ that threads the 
chain and couples to the fermionic fields in the form 
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c t C;e -<£^ 



Here A m m+ i is the vector 
potential defined on the link between sites m and m + 1, 
such that J2m=i A m ,m+i = The Hamiltonian in the 
presence of the electromagnetic field becomes 



H = E J * 

3<l 



1 



l\ (-l)'-j 



x 



(3.3) 



The current operator associated with a given link is de- 
fined as 



dH 



dAj 

££(-l)™J(n) x 



n-l 



(3.4) 



n>0 £=0 

xc t g»Ep='l" _1 (T«p-- A P,P+l) ( 



+h.c. 



Therefore the integrated current operator in the limit of 
weak fields, A/j+i — > 0, reads 



^ = E^ 



(3.5) 



^^(-l)"nJ(n)cj e »S 



j+n-l 
P=J 



r 



j n^O 



Reverting to spin operators, we obtain the spin current 
operator 



(3.6) 



Notice that J is a nonlocal operator, since it acts on sites 
separated by arbitrary distances. 

We can now compute the commutator of the current 
operator with the Hamiltonian in Eq. (3.1). The result 
is 



[J,H] 



- E t( m _ n )Jmn{Jln ~ Jim) 

+(2l-m~n)J lm J ln ]S[S+S~. (3.7) 



In the thermodynamic limit, J(n) — > 1/n 2 and we find 

[J,H}=0. (3.8) 

Therefore, the current operator is conserved and trans- 
port is ballistic in the thermodynamic limit. That the 
conservation law does not hold for a finite chain can be 
easily seen by considering a three site ring. In this case 
the Hamiltonian and current operator of the Haldane- 
Shastry model reduce to those of the Heisenberg model, 
and these two operators clearly do not commute. 



The current operator is orthogonal to the set of con- 
served quantities associated with integrability. For in- 
stance, the first nontrivial conserved quantity is 



Q 3 



; E 

ijk 



Zi Zj Zfz 
%ij Zjk%ki 



S t ■ (Sj x S fc ), 



(3.9) 



where Zj = e l27rj//L and Zij = z { — zj. The orthogonal- 



ity can be shown by arguing that the current operator 
is odd under a ir rotation about the x axis, which takes 
Sj —¥ —Sj and Sj — > — Sj , whereas the conserved quan- 
tities derived in Ref. [74] are SU(2) invariants and there- 
fore even under this transformation. Unlike the current 
operator, these conserved quantities commute with the 
Hamiltonian for a finite chain with periodic boundary 
conditions. 



IV. THE ATTRACTIVE HUBBARD MODEL 

We now want to consider the negative U Hubbard 
model at 1/2-filling 



H 



KiC aj +i + h.c.) - U{n 3 - 1/2)' 



(4.1) 



We may take the continuum limit and bosonize the 
resulting Dirac fermions, which now carry a spin index, 
finally introducing charge and spin bosons. These are de- 
coupled up to irrelevant operators. For U < and zero 
magnetic field, the spin excitations are gapped and we ex- 
pect the associated contribution to the low frequency con- 
ductivity to be Boltzmann suppressed at temperatures 
small compared to the spin gap. The low energy effec- 
tive Hamiltonian contains only the charge boson, which 
we again call (f> for simplicity. This low-energy Hamilto- 
nian now contains the Umklapp scattering term 



dx cos (v87r <?!>). 



(4.2) 



The bare coupling constant A is oc — U at small U. It is 
marginally irrelevant, as in the spinless model at A = 1, 
and the effective coupling at scale T behaves similarly to 
Eq. (2.24) 



1 lnA(T) , 
H = In 



X(T) 



T 



(4.3) 



where the cut-off scale T depends on U. Thus the self- 
energy for the charge boson has the form of Eq. (2.20), 
with 



2 7 (T) = ir\ 2 (T)T 



(4.4) 



as in Eq. (2.23). 

This lowest order, RG improved, calculation again ig- 
nores the possibility of a Drude weight. A rigorous Mazur 
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bound on the Drude weight for the Hubbard model was 
established in [7], using the conserved energy current, 
but it fails at half-filling where (JJe) = 0. In this case, 
other local conserved charges have not been considered. 

In the limit U/t — > — oo, the spin gap becomes infinite 
and the model is equivalent to the Heisenberg model with 
J oc t 2 /U and charge operators replacing spin operators. 
In this limit, we may import all results discussed above 
for the Heisenberg model. The fact that the Hubbard 
interaction is marginal (K — 1), as well as the equiva- 
lence with the Heisenberg model at large U follow from 
the exact duality transformation 



(4.5) 



which changes the sign of U. This maps bilinear opera- 
tors as follows 



TV 

h.c. -> 



-(c]a z c j+1 - h.c). (4.6) 



The charge current maps into the spin current. Adding a 
chemical potential to dope away from 1/2-filling at nega- 
tive U is equivalent to adding a magnetic field at positive 
U. In the large negative U limit, where we only allow 
empty sites and double occupancy, this becomes equiva- 
lent to the Heisenberg model (A = 1). The spin SU(2) 
symmetry of the U > model implies a "hidden" SU(2) 
symmetry in the charge sector. Thus, the Hubbard model 
at half-filling actually has SU{2) x SU{2)/Z 2 = S*0(4) 
symmetry for either sign of U . 

The continuum limit model also has a duality symme- 
try 



(4.7) 



Umklapp at U < maps into the ± component of the 
marginal spin operator at U > 0: 

^Ll^U^Rt^Rl -> ^Li^L^R^Ri- (4.8) 

Similarly, the product of left and right charge currents 
maps into the product of left and right z-components 
of spin currents: JlJr — > ^Jl^r- There is an exact 
SU(2) x SU (2) symmetry of the continuum model. Umk- 
lapp and JlJr interactions are sitting exactly on the sep- 
aratrix of the Kosterlitz-Thouless RG flow. 

We might think of a negative U as arising from phonon 
exchange. More general models, with longer range inter- 
actions will not be on the separatrix but may have the 
low energy Hamiltonian, with Umklapp and JlJr in- 
teractions. For a range of interaction parameters, the 
Umklapp will be irrelevant and we get a Lorentzian con- 
ductivity. 

The negative U Hubbard model may be thought of as 
being in a one-dimensional version of a superconducting 



state. However, it is important to realize that it is not 
a true superconductor and the presence or absence of 
a finite Drude weight at finite T is independent of the 
absence of true superconductivity. One way of seeing this 
point is to consider a ring of circumference L penetrated 
by a flux $. If the ring was a torus of sufficiently large 
thickness, we might be able to model its properties using 
London theory. Using the London equation 



n s {e*f 



mc 



(4.9) 



and A — <&/£, we obtain a persistent current ocl/L 



J = -- 



mcL 



(4.10) 



Here e* = 2e is the charge of the Cooper pairs. In London 
theory, this formula is true at finite T with n s (T), the 
superfluid density, a decreasing function of T which is 
non-zero for T < T c . On the other hand, for the one- 
dimensional Hubbard model, the finite size spectrum for 
current-carrying states is: 



7TO, 



E= ^p(n + $/27r) 2 + £ 



(4.11) 



where v c is the velocity of charge excitations. (This result 
is obtained ignoring the irrelevant Umklapp interactions.) 
The partition function is 



E 



exp[-(7ra c /£T)(n + $/2^) 2 ], (4.12) 



and the resulting persistent current 
j = dF/d<5>. 



(4.13) 



For sufficiently large L such that T is much greater than 
the finite size gap, T >• v c /L, this gives 



j ss — 2sin$ 



-kLT/v c 



(4.14) 



This is exponentially small in L, unlike the case for a 
superconductor, Eq. (4.10). The superfluid density van- 
ishes at any finite T. On the other hand, if we ig- 
nore Umklapp scattering, the model has a T-independent 
Drude weight D — v c /4ir, as in Eq. (2.31). 

>From the London equation for a superconductor it 
follows that 



lim to Im <t(uj 7 q) 



n s (e*f 



(4.15) 



independent of q. On the other hand, for the Hubbard 
model, Im a(q,uj) can be obtained from Eq. (2.32) with 
v — > v c and it is easy to see that lo Im a(q, ui) vanishes in 
the limit lo — > at non-zero q. 
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V. SUMMARY AND CONCLUSIONS 

To summarize, we have studied how nontrivial con- 
servation laws affect the transport properties of one- 
dimensional quantum systems. We focused, in particular, 
on the spin current in the integrable XX Z model. Away 
from half-filling, a part of the spin current is protected 
by conservation laws leading to a ballistic channel which 
coexists with a diffusive channel. For the half-filled case, 
however, none of the infinitely many local conservation 
laws responsible for integrability has any overlap with the 
current by symmetry. Nonetheless, several works have 
argued in favor of ballistic transport at finite tempera- 
tures even in this situation. This would require the ex- 
istence of an unknown nonlocal conservation law which 
has finite overlap with the current operator. To inves- 
tigate this controversial problem, our strategy was the 
following: Assuming that such a conservation law does 
not exist, we calculated the current relaxation within the 
Kubo formalism using a self-energy approach for the bo- 
son propagator. Since the parameters in the field theory 
are explicitly known due to the integrability of the micro- 
scopic model we obtained a parameter-free formula for 
the optical conductivity. Any shift of weight from this 
regular into a Drude part can then be described within 
the memory-matrix formalism and explicitly tested for 
by comparing the parameter-free result with numerical 
and experimental data. 

By a numerical study of the time-dependent current- 
current correlation function we have shown that the in- 
termediate time decay is well described by the calculated 
relaxation rate. This excludes, in particular, the large 
Drude weight found by Bethe ansatz in Ref. [9]. As we 
have shown, this result corresponds to our field theory 
with the relaxation rate set exactly to zero by hand. A 
small Drude weight can, however, not be excluded by 
the intermediate time data. This includes, in particu- 
lar, the Drude weight found in a different Bethe ansatz 
calculation. 8 At least at high temperatures it is, how- 
ever, known that this approach violates exact relations 
so that the solution cannot be exact. It is not clear at 
present if or why these results could be considered as 
an approximate solution at low temperatures. Quantum 
Monte Carlo calculations, on the other hand, are per- 
formed using imaginary times. The results presented in 
Refs. [11,12] therefore cannot be used to resolve a decay 
rate smaller than the separation of Matsubara frequen- 
cies. An interesting possibility is the idea to analyti- 
cally continue our results to imaginary frequencies and 
to check if this is consistent with Quantum Monte Carlo 
data. Such an analysis has been performed recently for 
the isotropic case and good agreement was found. 30 It 
would be very interesting to perform a similar analysis in 
the anisotropic case for A values such that the decay rate 
is not too small. The advantage would be that additional 
complications due to logarithmic corrections present at 
the isotropic point do not occur. In Ref. [19] we also pre- 
sented numerical data for the current-current correlation 



funciton at infinite temperatures and showed that even 
in this large time scale persists. It is therefore 

unclear how data obtained by exact diagonalization can 
be reliably extrapolated to the thermodynamic limit. 

>From our perturbative result for the boson prop- 
agator at finite temperatures we found that the long- 
wavelength contribution to the spin-lattice relaxation 
rate is diffusive. By comparing with experiments on the 
spin chain compound S^CuOs we showed that our for- 
mula describes the experimental data well. Importantly, 
we showed that the data for the magnetic susceptibility 
and the q ~ and q ~ tt contributions to the spin-lattice 
relaxation rate are all consistently described by the field 
theory formulas with the same exchange constant J and 
using the hyperfine coupling tensor as determined exper- 
imentally. We also pointed out that our results are con- 
sistent with a previous theory of electron spin resonance 
in spin chains. 60 For the longitudinal spin structure fac- 
tor at zero magnetic field our theory predicts a crossover 
from the known non-Lorentzian lineshape with linewidth 
~ q 3 at zero temperature, to a Lorentzian lineshape with 
a linewidth set by the relaxation rate 7 at sufficiently 
high temperatures. 

We also discussed the integrable Haldane-Shastry 
model where the spectrum is known to consist of free 
spinons in contrast to the interacting spinons in the 
Heisenberg model. As in the Heisenberg model at zero 
magnetic field the current operator does not have any 
overlap with the known conserved quantities. However, 
in this case the nonlocal current operator itself becomes 
conserved in the thermodynamic limit and transport is 
therefore ballistic. 

Finally, we showed that our calculations for the spin 
current in the XX Z model carry over to the charge cur- 
rent in the attractive Hubbard model at half-filling. We 
stressed the point that an infinite dc conductivity does 
not imply that the system is a true superconductor. This 
can be seen by calculating the superfluid density which 
turns out to be zero in this case. 



Appendix A: Spin Green's function 

Consider a boson propagator with a momentum- 
independent decay term (see Eq. (2.16)) 

poo poo 

Xret(q,u) = -if dt dx 

Jo J -00 

x e j ^-^)([a^(t,x),a^(o,o)]) T 



L J 2 -q 2 + 2l 1 (T)u} , y ' 

for small oj and q. (We set the velocity equal to 1. 7 must 
be > by causality.) We wish to calculate the long-time 
behavior of the correlation function 

G(x,t) = ^-(d x <j>(t,x)d x <j>(0,0)) T , (A2) 
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which can be expressed in terms of Xret{l,u) with the 
help of Eq. (2.56). Using this equation and (Al) we find 



■ [ diUdq e -i(ut-qx) q 2 

{X ' t} ~* ~ l J (2ir) 2 (1 - e-l^)[q - (cu 2 + 2i 7 u)i/2] [q + (lu 2 + 2i 7 o,) */2] ~ ^ ^ ^ ^3) 

We now consider only the first term; we return to the 7 — > —7 term later. We are defining (lu 2 + 2i^Lu) 1 l 2 to be a 
particular branch of the square root; let us specify carefully which one. It will be convenient to define this branch 
with the branch cut along the negative imaginary oj axis between lu — and lu = —2i-f. Thus 

(lu 2 + 2ijLu) 1/2 -> - u 2 (uj = -iu + 5, < u < 2 7 , (5^0+) 

-> lo + i"f (\lu\ > 7) 

-> -y/2ju - u 2 (w = -iu-<J, < u < 27, <5^0 + ). (A4) 



We do the g-integral first. Noting that the result is man- 
ifestly an even function of x, we just consider explicitly 
the case x > 0. Noting that this branch of [lj 2 + 2i^Lu) 1 / 2 
always has a positive imaginary part for real lu, we may 
close the q integral in the upper half plane, encircling the 
pole at q = {lo 2 + 2i^lu) 1 / 2 , giving 

— exp[-icut + i(uj 2 + 2iju) 1/2 \x\] 

Consider the analytic structure of the integrand in the 
complex lu plane. (We consider only the first term; we 
return to the 7 — >• —7 term later.) There is a branch 
cut along the negative imaginary axis from lu = to 
lu = —2i-f. There are also poles at uj = 2mnT for n = ±1, 
±2, . . .. The integrand behaves as l/cu 1 / 2 at lu — > so the 
integral is convergent. Let's assume \t\ > \x\. Then, we 
can close the integral in the upper half plane for t < or 
the lower half-plane for t > 0. First consider the simpler 
case t < 0. Then we only pick up the contributions from 
the poles lu = 2mnT, n = 1, 2, 3, . . . Let us now assume 
7(T) < T. Note that we expect this to be true at low T, 
7(T) oc T 4K - 3 , K > 1. Then, at the poles: 



(lu 2 + 2i~/Lu) 1/2 w lu = 2mnT. 



(A6) 



Then, for t < 0, the first term gives the approximately 
7-independent result: 

00 

G first *) ^ _7rT2 S n e^^-Vi+W^T)!,!) (A?) 

71=1 

_ e -7kl 1 



4tt {(l/7rT)sinh[7rr(t - \x\ 



12 • 



Actually, while we can always ignore the 7-dependent 
corrections to the pre-factor, those corrections in the ex- 
ponent cannot be ignored at sufficiently large \x\. We 



have Taylor expanded the square root in the exponential 
to first order in 7/T. We see that the correction becomes 
important at \x\ of order I/7. At such large values of |x|, 
higher order corrections must also be included. If we 
assume \x\ <C I/7 we may simply drop the e -7 '^' factor. 

Now consider the first term in Eq. (A5) for t > 0. 
There are now both pole and cut contributions. For the 
pole contributions we find again 



Qt>0 

first, pole 



- 7 |x| 



4tt {(1/ttT) sinh[7rT(t - 



(A8) 

But now, we must also consider the cut contribution. 
We may Taylor expand the denominator in Eq. (A5) to 
first order since 0\cj\ <C 1 along the cut, lu = —iu, with 
< u < 2 7 < T. Thus: 



Qt>0 

first, cut 



> -y / 2 1 u-u 2 e- ut 



(A9) 



1 /•-' <1 

47T/3 J u 

x{exp[— i\p2rfu — u 2 \x\\ + cxp[i\/2ju — u 2 |x|]}. 



Now let us assume \t\ ^> I/7, so that the integral is 
dominated by u <C I/7 and we may approximate it by: 



Qt>0 

first, cut 



x cos 




(A10) 



Here we have expanded the square root to first order 
in the exponential. This first correction, as well as the 
higher order ones, may be dropped provided that: 

N<|t|V7|tf. (All) 

Note that this is automatically true since we have already 
assumed |x| < \t\ and y/~f\t\ 3> 1. Assuming Eq. (All) 
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and letting u — v 2 /2, this becomes Now consider the 7 — > —7 term in Eq. (A3). We define 

^ a convenient branch of \Juj 2 — 2iju> with the branch cut 

(jt>o _j. /" rf-y e _t,2 l*l' /2 [e _i ^ 7 ' :l: ' t ' + c c 1 along the positive imaginary axis from u) — to us = 2i-f 

first,cut 2tt/3 7 L ' ' J obeying: 



= T,fe e -^ 2 /( 2 l*D. (A12) 



(uj 2 -2i-fcj) 1/2 ->■ V27M - u 2 (w = m + <5, 0<w<27, <5 ^ 0+) 
->• w - «7 (|w| > 7) 

-> -V/27U-M 2 (w = m-(5, 0<u<2 7 , (5^0+) (A13) 



Now, choosing again a; > 0, the (7 integral, in the upper 
half plane, encircles the pole at q = —(lu 2 — 2ijuj) 1 ^ 2 , 
giving: 

CsecondM) -> / ^ exp[-iwi - i(uj 2 - 2i 7 w) 1 / 2 |x|] 
Now the pole terms give 



G a 



*second,pole 4 ^ {(i/^ T ) sinlnVr(t + |x|)]} 2 

(for either sign of t). The cut term, is now present for 
t < 0, giving: 



G *<o =G t>0 > Ti / 27r7 c ^7^ 2 /(2|t|) 

second, cut first, cut \/ k| 



(A16) 



Combining both terms gives 

GOM) ->G (a:,i) + G int (a;,t), (1*1 > M)- (A17) 
where Gq(x, t) is the result in the non-interacting case 



G (x,t) 



+ 



1 f e-Tl^l 
4^r \{(1/ttT) sinh[7rT(t-x)]} 2 

e 7k| 

{(1/ttT) sinh[7rT(f + a;)]} 2 



(A18) 



and 



G int (t,x)^TJ 2 -^e-^ 2 ^ (A19) 



(for both signs of t). 

Note that we assumed \t\ > \x\ below Eq. (A5) deter- 
mining in which half-plane the contour for the w-integral 
was closed. We may attempt to evaluate the integrals 
using the same method also when |a;| > \t\. Now the u> 
integral is closed in the upper half plane for the first term 
and the lower half plane for the second, regardless of the 
sign of t. Thus the cut terms do not appear and we only 
get the pole terms: 



G(x,t) -> G Q (x,t) (\t\ < \x\). 



(A20) 



Let's reiterate the various conditions on t and x in 
order for Eqs. (A17) and (A20) to hold. The form of G 
depends on the assumption \x\ <C I/7 but not an any 
particular assumption on t. On the other hand, the form 
of Gint, which is only present for \t\ > \x\, also assumes 
1*1 > 1/7- 
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